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ABSTRACT 

Context. Previous far-infrared observations at low-angular resolution have reported the presence of water associated 

with low-velocity outflow shocks and protostellar envelopes. 

Aims. We want to elucidate the origin of water emission in protostellar systems. 

Methods. The outflow driven by the intermediate-mass class protostar CepE is among the most luminous outflows 
detected so far. Using the IRAM 30m telescope, we searched for and detected the P-H2O 3i3 — 220 line emission at 
183 GHz in the CepE star-forming core. The emission arises from high-velocity gas close to the protostar, which is 
unresolved in the main beam of the telescope. Complementary observations at 2" resolution with the Plateau de Bure 
interferometer helped establish the origin of the emission detected and the physical conditions in the emitting gas. The 
water line profile and its spatial distribution are very similar to those of SiO. We find that the H2O emission arises from 
warm (~ 200 K), dense ((1 — 2) x 10® cm~'^) gas, and its abundance is enhanced by one to two orders of magnitude 
with respect to the protostellar envelope. 

Results. We detect water emission in strong shocks from the high-velocity jet at 1000 AU from the protostar. Despite 
the large beam size of the telescope, such emission should be detectable with Herschel. 
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1. Introduction 

Water is a key molecule not only for oxygen chemistry in 
the interstellar gas, but also in the dynamical evolution 
of star- forming regions, because it is one of the main gas 
coolants (Kaufman & Neufeld, 1996; Bergin et al. 1998). 
Whereas its abundance is very low in dark clouds, prestel- 
lar cores, and starless regions (from less than 10~^ to a few 
lO"'^; Bergin & Snell, 2002; Moneti et al. 2001), it can be 
greatly enhanced in star-forming regions, through the sput- 
tering of frozen water from grain mantles and through high- 
temperature, sensitive reactions in the gas phase (Elitzur & 
de Jong, 1978; Ceccarelli et al. 1996; Kaufmann & Neufeld 
1996). It has been proposed that such processes are driven 
by the heating of the nascent stars (hot core/corino re- 
gion) or by shocks either in the accretion region in the in- 
ner protostellar core or in the outflow carrying away the 
angular momentum of the envelope. Previous observations 
of water lines obtained with the Infrared Space Observatory 
have substantially confirmed these theoretical expectations 
(see e.g. Giannini et al. 2001; Ceccarelli et al. 2000; van 
Dishoeck 2004 for a review). However, many details are 
still missing and the relative efficiency of these mechanisms 
is not known, so that the origin of the water emission in 
protostellar regions remains a puzzle to be solved. One ob- 
stacle is indeed the relatively low spatial resolution of these 
observations. 

Most of the H2O lines cannot be observed from ground 
because of strong atmospheric absorption. The para-H20 
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line 3i3 — 220 at 183.3 GHz is an exception since it can 
be enhanced by the masing effect even in regions of moder- 
ate excitations. Using the IRAM 30m telescope, Cernicharo 
et al. (1990,1994) detected large-scale P-H2O 3i3 - 220 hne 
emission in giant molecular clouds and star-forming regions 
and they showed that this line can be used to investi- 
gate the spatial distribution of H2O in protostellar envi- 
ronments at an angular resolution comparable to the size 
of the protostellar cores (14" at the IRAM 30m telescope). 
Observing the low-mass protostellar outflows HH 7-11 and 
L1448, Cernicharo et al. (1996) detected water line emis- 
sion arising from shocks associated with the Herbig-Haro 
objects in the outflows. 

In this Letter, we report on observations of the 
P-H2O 3i3 — 220 hue in the star-forming core CepE with 
the IRAM 30m telescop^ Located at 730pc, it is a rela- 
tively nearby intermediate-mass Class source of 100 Lq 
and a mass of 18 Mq, intensively studied at mm and IR 
wavelengths (Moro-Martin et al. 2001; Noriega-Crespo et 
al. 2004; Lefloch et al. 1996, hereafter L96). The protostar 
drives an exceptionally powerful outflow, whose southern 
lobe is terminated by the Herbig-Haro object HH 377. Far- 
infrared observations at low spectral and angular resolu- 
tion with ISO have revealed large amounts of water, which 
Giannini et al. (2001) interpret as arising from low-velocity 
shocks, whose far-infrared luminosity amounts to L ~ 3 Lq. 



^ Based on observations carried out with the IRAM 30m 
Telescope. IRAM is supported by INSU/CNRS (France), MPG 
(Germany), and IGN (Spain) 
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Fig. 1. (left) Map of the P-H2O 813 — 220 emission between 
— 110 and +70 km s^^ in the outflow of Cep E obtained with 
the IRAM 30m telescope. Coordinates are in arcsec off- 
sets with respect to the nominal position of the protostar. 
(right) H2O and SiO J = 8 — >■ 7 emission detected towards 
the protostellar source. The H2O emission is averaged over 
a region of 6" x 6" around the protostar. Intensities are 
expressed in units of antenna temperature. The different 
components of the spectrum are marked with arrows. 



Observations of the millimeter lines of SiO observa- 
tions reveal the presence of strong shocks along the outflow 
(L96). Because it is usually undetected in the cold, qui- 
escent molecular gas, SiO is a particularly good tracer of 
shocks that arc strong enough to release refractory elements 
from grain cores. Recent Herschel observations suggest that 
SiO is also a good tracer of the water-emitting shock regions 
(Lefloch et al. 2010; Nisini et al. 2010). We present obser- 
vations of the high-excitation SiO J = 8 — > 7 line towards 
the protostar. 
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Fig. 2. (left) Map of the velocity-integrated SiO emis- 
sion in the Cop E jet obtained at 2" resolution with the 
Plateau de Bure interferometer (PdBI; Lefloch et al. in 
prep.). The SiO J — 2 1 emission integrated between 
-1-40 and -f-80kms~^ (—90 and — 40kms~^) is drawn in 
solid (dashed) contours, superposed on the 3 mm con- 
tinuum emission map (grayscale and white contours). 
First contour and contour interval of the jet emission are 
10-2 and 5 x lO'^ (lO'^) Jy/beam for the redshifted 
(blueshifted) emission. First contour and contour inter- 
val of the 3mm continuum emission is 2 x 10~"^ Jy/beam. 
The main beam of the IRAM 30m telescope; at the fre- 
quency of the P-H2O 3i3 — 220 line is marked by a circle 
(dashed). The synthesized beam of the PdBI (2.4" x 2.0") 
is indicated as a filled ellipse in the bottom right cor- 
ner, (right) Comparison of the PdBI SiO J = 2 —> 1 and 
the average P-H2O 813 — 220 spectra towards the protostar. 
Intensities are expressed in units of antenna temperature. 



2. Observations 

2.1. Water observation 

We observed the line of para-H20 3i3 220 at 183.3 GHz 
at the IRAM 30m telescope in January 2009, near the tran- 
sit of CcpE between 61° and 65° of elevation. Both receivers 
at 2mm were used simultaneously. Observations were car- 
ried out in wobbler switching mode with a throw of 180". 
Filterbanks with IMHz spectral resolution and an autocor- 
rclator providing a spectral resolution of 1.25 MHz (2 km/s) 
were used as spectrometers. Weather conditions were very 
good and stable, with an atmospheric optical depth « 1.10 
along the line of sight. Calibration was carefully monitored 
every 4 minutes between two integration scans on CepE. 
Water vapor and system temperatures were found to re- 
main stable, in the range 1650 — 1800 K. Any fluctuation of 
the sky emissivity between the on and off positions would 
result in a poor spectroscopic baseline and detection of the 
mesospheric water line, which is not the case in the present 
observations. Pointing was checked against nearby quasars 
every hour and has an accuracy of « 2". At the frequency 
of the H2O line, the main-beam efficiency of the telescope 
is 0.63 and the half-power beamwidth is 13.2". We did a 
small oversampled map at 6" sampling of the southern lobe 
of the outflow, centered on the nominal position of the pro- 



tostar 23'*03"13.0" +61°42'21" (J2000). The typical rms is 
lOOmK per velocity interval of 2 km s~^. 



2.2. SiO observations 

The SiO J = 8 -J> 7 transition at 347.33006 GHz was ob- 
served towards the protostar at the; JCMT in August 2005. 
The weather conditions were good with an average sys- 
tem temperature Tsys ~ 400 K. Observations were carried 
out with a nutating secondary, adopting a throw of 60". 
The DBS was used as spectrometer, providing a spectral 
resolution of 625kHz. The resolution was degraded to 3.75 
MHz, corresponding to a flnal velocity resolution of about 
3.2 km s~^. At the frequency of the SiO line, the main-beam 
efficiency of the telescope is about 0.60 and the half-power 
beamwidth is ~ 14.2", very similar to the telescope param- 
eters of the P-H2O observations. Complementary interfer- 
ometric observations at 2" resolution with the Plateau de 
Bure in the SiO J = 2 — >^ 1 transition provide us with the 
spatial distribution of the emission and its excitation con- 
ditions. These observations will be presented and discussed 
in detail in a forthcoming paper (Lefloch et al. in prep.). 
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3. Results and discussion 

3.1. Water emission in CepE 

The map of the P-H2O 813 — 220 emission in CepE is dis- 
played in the left panel of Fig. 1. Emission is detected only 
towards the protostellar core, in a region unresolved by the 
main beam of the IRAM 30m telescope. In particular, we do 
not detect any H2O emission towards HH77, the southern 
bowshock at offset position (—5", —15"), associated with 
bright emission in the lines of SiO and H2 (L96). To im- 
prove the signal-to-noise ratio of the data, we averaged the 
emission of the water line over the inner core, in the range 
Aa = [1"; -5"] A6 = [1"; -5"]. The averaged spectrum is 
displayed in the right panel of Fig. 1, and the final rms 
noise is 30mK (TJ) per interval of 2 kms"-'^. 

Two bright components 1 K) with broad linewidths 
(~ lOkms^^) are detected at about visr ~ —65 
and — 45kms~^ in the direction of the protostar. Two 
other, weaker, H2O components are detected at approxi- 
mately symmetric (redshifted) velocities visr ~ +34 and 
-|-50kms~^. We do not detect any emission at the veloc- 
ity of the protostellar envelope (u/sr = — 10.9kms~^; L96), 
hence excluding any contribution from the hot corino to the 
water emission. The presence of both blue- and redshifted 
components in the water line spectrum indicates that the 
emission arises from the central protostellar region. This is 
because there is barely any overlap in the plane of the sky 
between the blue and redshifted wings of the Cep E outflow, 
except in the central protostellar region as a consequence 
of the limited angular resolution of the observations (see 
Figs. 1-2; also L96). 

The spectrum of the SiO J = 8 ^ 7 emission in the di- 
rection of the protostar is displayed in Fig. 1 and looks very 
similar to the P-H2O 813 — 220 line spectrum. Resemblance 
between both emissions is even more striking when com- 
paring the water spectrum with the J = 2 — 1 emission 
observed with the PdBI (Fig. 2). This suggests that the 
emission of both tracers arises from the same region. 
Observations of the J = 2 — ^ 1 transition with the PdBI of 
the central protostellar regions (Lefloch et al. in prep) show 
that the SiO emission arises from two compact clumps, 
in the blue- and redshifted gas, respectively, located at 
« 1.5" from the protostar. They have a typical size of 2.5" 
(2.3 X 10^^ cm), and a transverse size of 1" (Fig. 2). Their 
similar physical properties (size, velocity, distance) suggest 
they are related to an episodic ejection phenomenon. In 
what follows, we assume that both the II2O and the SiO 
emissions have the same size. 

3.2. Physical conditions 

The gas kinetic temperature Tk in the redshifted, high- 
velocity jet was previously determined by Hatchell et al. 
(1999) from JCMT observations of the CO J = 4 ^- 3 and 
J = 2 — > 1 lines. Unfortunately, their analysis is biased by 
a highly overestimated value of the main-beam efficiency 
TiMB of the J = 4 — 3 observations, poorly known at that 
time, leading to an underestimate of Tk (~ 50 K). Based 
on their observations and adopting a standard value con- 
sistent with the subsequent determinations from point-like 
calibratorf]^ '7ms — 0.35, we obtain the jet kinetic temper- 
ature and column density, Tk — 120 K and N(CO) ~ 2.5 x 

^ http:/ /www.jach. hawaii.edu/JCMT/spectral_line/Standards/ 



10^^ cm^^, respectively (see also L96). Adopting a stan- 
dard CO abundance in the gas, we then derive the source- 
averaged gas column-density N(H2) ~ 2.5 x lO'^^ cm^^. 
We adopt these values of Tk and N(CO) for both water- 
emitting clumps at blue- and redshifted velocities. 

The gas density and column density in the water- 
emitting clumps were derived from modeling the II2O and 
SiO line brightness temperatures in the large-velocity gra- 
dient (LVG) approach. We adopted the collisional coeffi- 
cients of Faure et al. (2007) and Dayou & Balanga (2006), 
respectively. The gas density can be constrained reason- 
ably well from the ratio of the SiO J = 8— >-7toJ = 3— >-2 
brightness temperatures. This ratio was computed using 
the SiO J = 3 — >■ 2 observations in the Cep E protostellar 
core obtained by L96 (see also Fig. 3). We note that the 
source is unresolved by the single-dish JCMT and IRAM 
30m observations; therefore, we approximated the ratio of 
the line brightness temperatures by the ratio of the main- 
beam brightness temperatures corrected for the difference 
of main-beam solid angles R(8 — 7/3 — 2). A description 
of the radiative transfer code and the procedure applied to 
deal with the masing effect is presented in Cernicharo et al. 
(1994) and Gonzalez-Alfonso & Cernicharo (1993). 

The top right panel of Fig. 3 displays the variations of 
R(8 — 7/3 — 2). The ratio varies little in the low- velocity 
gas and is about 0.2. It increases as a function of velocity 
and reaches values up to 0.7 and 1.5 in the redshifted and 
blueshifted gas, respectively. The variations of R(8— 7/3— 2) 
as a function of T^ and n(H2) were computed for a typical 
column density N(SiO) = 2 x lO^** cm~^ and a linewidth 
Ay — lOkms^^, and they are displayed in Fig. 3. With 
the constraint Tk > 100 K, it comes out that n(Il2) must 
be about 1.0 x lO'' cm"'^ in order to account for the line 
intensity ratio in the redshifted gas . The density in the 
blueshifted gas is slightly higher « 2 x 10® cm~'^. These cal- 
culations show that the ratio depends weakly on Tk above 
100 K. For these values of n(H2) and the kinetic tempera- 
ture estimated above, we estimate the source-averaged col- 
umn density necessary to account for the J = 8 — 7 line 
intensity N(SiO) = (2.0 - 3.3) x W'^ cm'^ for gas at 120 K 
in the blue- and redshifted components. Under such phys- 
ical conditions both SiO transitions are optically thin, so 
that R(8-7/3-2) depends weakly on the adopted value of 
the gas column density. From comparison with N(CO) and 
assuming a standard abundance [C0]/[H2] = 10"'', we ob- 
tain the relative SiO abundance in the shock : [SiO]/[Il2] ~ 
1.0 X 10"^, similar to the abundances measured in other 
Class protostellar outflows (Lefloch et al. 1998). 

The para- water (source-averaged) column density N(p— 
H2O) was determined from modeling the line brightness 
temperature Tb- The latter was derived from the main- 
beam brightness temperature T^v/b = Tb x //, where the 
filling factor // was obtained assuming that the H2O and 
SiO J = 2 — ^ 1 emitting regions have similar sizes, in agree- 
ment with Sect. 3.1. We find Tb = 60-80 K and 15-29 K in 
the blueshifted and redshifted gas, respectively. The bright- 
ness temperature of the P-H2O 813 — 220 line was computed 
for a wide range of physical conditions, adopting a typi- 
cal linewidth Aw = lOkms"^ : Tk between 50 and 200 K, 
n(H2) in the range 10^ - 10^ cm"^ and N(p - H2O) in the 
range 3 x 10^® to 10^^ cm~^. The results of the calculations 
are presented in Fig. 4. 

It comes out that the kinetic temperature Tk has to 
be higher than 100 K to account for brightness temper- 
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Fig. 4. Predictions of the P-H2O 813 — 220 line intensity, in 
units of brightness temperature, for gas with density n(H2) 
in the range 10^ — 10* cm^"^, temperature in the range 50 — 
200 K, for column densities of para-water equal to 0.003 
(bottom curve), 0.01, 0.03, 0.1, 0.3, 1.0 x lO^^ cm'^ (top 
curve). The range of Tb — n(H2) determined in the water- 
emitting clumps is delineated by rectangles. 



atures of 60 — 80 K. Our calculations indicate a P-H2O 
column density in the range (0.3 — 2.0) x 10^^ cm~^ for 
Tk — 200 K. Assuming an ortho-to-para ratio of 3, we de- 
rive N(H20) = (1.2 - 8.0) X 10" cm-2, and using the H2 
column density derived in the high-velocity clumps, we ob- 
tain the water abundance X(H20) = (0.5 — 3.2) x 10^*. In 
this range of physical conditions, the line optical depth r is 
small w — 1 and shows evidence for the weak masing ampli- 
fication of P-H2O 3i3 — 220 line. Suprathermal excitation of 
the line would require much higher H2O column densities. 
This is possible if the emission arises from a small region of 
subarsec size (see e.g. Cernicharo et al. 1994) and is very 
sharp. The line Tb would then look weaker as a result of 
the dilution in the 2kms~^ channels of the spectrometer. 

Moro-Martin et al. (2001) analyzed the far-infrared 
emission lines of CO and H2O in CepE, as observed with 
ISO. Like Giannini et al. (2001), they found that the line 
emission can be accounted for by an extended component of 
low-density gas (2x lO'* cm^^) at ~ 1000 K. They show that 
the CO and H2O lines could be reproduced equally well by 
"warm" gas at a much lower temperature (~ 200 K) and 
higher density (2 x 10^ cm~^). The physical conditions we 
derive for the P-H2O 3i3 — 220 emitting region are actually 



in rather good agreement with this "warm" solution and 
inconsistent with the low-density "hot" solution. 

3.3. Shocked water in CepE 

Crimier et al. (2010) determined the physical conditions 
in the protostellar envelope of Cep E envelope and found 
n(Il2) = 1.0 X 10''cm~^ and Tk ~ 50 K, respectively, at 
the location of the water-emitting clumps, i.e. 1000 AU 
(1.4") from the protostar. No P-H2O 813 — 220 emission is 
detected from the envelope at the 2a level = 3K (Tb). 
From Fig. 4, we conclude that A^(Il20) < 10^^ cm~^ in the 
inner protostellar envelope, and the water abundance is en- 
hanced by 1 to 2 orders of magnitude in the high-velocity 
clumps. The authors estimate a radius of « 220 AU{0.3") 
for the hot corino region, and the density is predicted very 
high ^ 10* cm~^; again, no masing effect is to be expected, 
in agreement with our observations. 

The water-emitting clumps are located along the high- 
velocity jet, which suggests that they are closely related 
with the jet (Fig. 2). One possibility is that these clumps 
are tracing the jet /envelope interaction. However, the den- 
sity of the protostellar envelope, as estimated by Crimier 
et al. (2010) at the location of the water emitting-clumps 
(n(Il2) = 1.0 X 10^ cm^-^), is about one order of magni- 
tude higher than the density measured in the shocked gas 
(10^ cm"-^), in apparent contradiction with this hypothe- 
sis. The structure of the envelope could be actually more 
complex than the simple power-law distribution derived 
by Crimier et al. (2010), as proposed that multiple proto- 
stars are embedded inside the core (Ladd & Hodapp, 1997). 
Direct observational evidence for several condensations is 
still missing, however. Alternatively, we speculate that the 
water emission could arise from internal shocks in the en- 
trained material of the flow. Interferometric observations at 
subarcsec scale in the inner core (r < 1000 AU) are needed 
to determine the gas density and velocity field in order to 
determine the location and the origin of the water-emitting 
shocks, as well as their relation with the jet. 

Simple calculations in the LVG approximation predict 
that the H2O emission from the high- velocity clumps should 
be easily detected in several ortho- and para-Il2 transi- 
tions with the high-resolution spectrometer HIFI onboard 
Herschel. Such observations will allow us to better constrain 
the shock physical conditions (see e.g. Flower & Pineau des 
Forets, 2010). 

Acknowledgements. We thank an anonymous referee and the Editor 
for comments that helped to improve the manuscript a lot. This 
project was supported by the research grant SAB2009-0011 of the 
Spanish Ministry of Education. 



References 

Bergin, E.A., Neufeld, D.A., Melnick, G.J., 1998, ApJ, 499, 777 
Bergin, E.A., Snell, R.L., 2002, ApJ, 581, L105 

CeccareUi C, HoUenbach D., Tielenas A.G.G.M., 1996, ApJ, 471, 400 
CeccareUi C, Castets A., Caux E., et al., 2000, ApJ, 355, 1129 
Cernicharo J., Thum C, Hein H., et al., A&A, 1990, 231, L15 
Cernicharo J., Gonzalez- Alfonso, E., Alcolea, J., 1994, 432, L59 
Cernicharo J., Bachiller R., Gonzalez-Alfonso, E., 1996, A&A, 305, 
L5 

Crimier N., CeccareUi C, Alonso-Albi T., et al., 2010, A&A, 516, 102 
Dayou. F., Balanga, C, 2006, A&A, 459, 297 
Eisloeffel, J., Smith, M.D., Davis, C.J., 1996, AJ, 112, 2086 
Elitzur, M., de Jong, T., 1978, A&A, 67, 323 

Faure, A., Crimier, N., CeccareUi, C, et al., 2007, A&A, 472, 1029 
Flower, D.R., Pineau des Forets, G., 2010, MNRAS, 406, 1745 
Giannini, T., Nisini B., Lorenzetti D., 2001, ApJ, 555, 40 
Kaufman, M., Neufeld, D., 1996, ApJ, 456, 611 



B. Lefloch et al.: Shocked water in the CepE protostellar outflow 




-100 -50 50 
V (km/s) 



5.5 , 6 ,„6.5 
log nlHgJ 



Fig. 3. (left) SiO emission observed in the J = 3 — 2 (bot- 
tom) and J = 8 — 7> 7 transitions towards the protostellar 
source. Fluxes are expressed in units of main-beam tem- 
perature, (right) Observed variations of the ratio of the 
CO J = 8— >7toJ = 3— )-2 brigthness temperatures R(8- 
7/3-2) as a function of velocity (top) and predicted vari- 
ations as a function of the gas density ri(H2) and the ki- 
netic temperature Tk (bottom) for a typical column density 
N(SiO) = 2 X 10^'' cm-2. The contour level R= 0.7, 1.5 are 
drawn with a solid black line; the contours R= 0.5, 1.8 are 
drawn with a dashed red line and account for an uncertainty 
of 30%. 
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